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Abstract 

We define new observables sensitive to BFKL dynamics in the con¬ 
text of multijet production at the large hadron collider (LHC). We pro¬ 
pose the study of the inclusive production of three jets well separated 
in rapidity from each other, with two of them being very forward. We 
show that the tagging of a third jet in the central region of rapidity 
allows for a very strong test of the BFKL formalism. In particular, 
we have studied two projections on azimuthal angles for the differ¬ 
ential cross section which allow for the definition of many different 
observables whose behavior when varying the pt and rapidity of the 
central jet is a distinct signal of BFKL dynamics. In order to reduce 
the theoretical uncertainties and influence of higher order corrections, 
we propose the study of ratios of correlation functions of products of 
cosines of azimuthal angle differences among the tagged jets. 


1 Introduction 

The large number of events already recorded and those to be produced 
in the near future at the Large Hadron Collider (LHC) offer a unique 
opportunity to disentangle the region of applicability of asymptotic 
calculations of scattering amplitudes in the high energy Regge limit. 
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In this work we focus on the description of new observables where 
the Balitsky-Fadin-Kuraev-Lipatov (BFKL) formalism, at leading (LL) 
[H El El m El |6] and next-to-leading (NLL) [7118] accuracy, should ap¬ 
ply. In a nutshell, when calculating a scattering amplitude within 
this approach we single out those contributions at each order in per¬ 
turbation theory with the largest numerical value when the center- 
of-mass energy, ^/s, is asymptotically higher than any of the other 
Mandelstam invariants. These enhanced contributions are linked to 
the rapidity dependence of the observable under consideration. From 
a phenomenological perspective, it is important to find the window of 
applicability for this formalism. This means to identify observables 
where the BFKL approach is distinct, i.e., quantities where it fits the 
measured data and all other possible approaches (fixed order or other 
resummations implemented in general Monte Carlo event generators) 
fail. 

So far, the search for BFKL effects has had the general drawback 
of having collisions with too low ^/s or rapidity differences among the 
tagged particles in the final state. A further problem has been to con¬ 
sider observables which are too inclusive as to be able to claim that the 
cross section under study could only be described by BFKL dynamics 
and nothing else. A canonical example is the growth of the hadron 
structure functions at low values of Bjorken x in Deep Inelastic Scat¬ 
tering. Indeed it is possible to get a good fit of the combined HERA 
data for F 2 ^l with a NLL BFKL calculation {e.g., PEU]). However, 
it is equally possible to ht these data with other approaches. We must 
find other, less inclusive, observables to test small x resummations. 

The LHC solves these problems since the available energies are 
much higher than at the Tevatron or HERA and there is enough 
statistics to allow for the study of very exclusive quantities, with strong 
kinematical cuts. The experimental challenges in this direction are re¬ 
markable since it is needed a large rapidity span in the hnal states and 
a good resolution in azimuthal angles. We should not only consider 
the usual “growth with energy” signal, associated to the exchange of 
a hard pomeron, but also other footprints related to energy flow and 
azimuthal angle dependences. The latter are the main subject of the 
work here presented. 

Typical BFKL observables at the LHC are the azimuthal angle 
((/>) decorrelation of two tagged forward jets widely separated in ra¬ 
pidity, Y, with associated inclusive mini-jet radiation, the so-called 
Mueller-Navelet jets m- This multiple emission appears as a fast 
decrease of (cos(n^)) as a function of Y [T2l [T3l [HI [15]. How¬ 
ever, these differential distributions suffer from a large influence of 
collinear regions in phase space imiiTj. This is due to the fact that 
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<cos(n(/>)>~ exp (asy(xn(l/ 2 ) — xo(l/ 2 ))), where Ug is the strong 
coupling and Xn{l) is, in Mellin space, the n-th Fourier component 
in (j) of the BFKL kernel where the region 7 ~ 1/2 dominates for 
large Y. The n = 0 component is very sensitive to collinear dynam¬ 
ics well beyond the original multi-Regge kinematics. Even though 
it is possible to resum these collinear contributions “on top” of the 
BFKL original calculation we believe that it is more important at 
present to fix the real region of applicability of the original BFKL 
formalism by using observables which are far less sensitive to this 
collinear “contamination”. It is our target to find in this way distinct 
BFKL observables. An important step in this direction was taken 
in Ha IT] where it was proposed to remove the n = 0 dependence 
by studying the ratios Cm,n = {cos {m c^))/{cos {n (p)) which behave 
like ~ exp (a 5 K(xm(l/ 2 ) — Xn(l/2))). It is important to note that 
the BFKL kernel for n 7 ^ 0 is insensitive to collinear regions, as it 
was shown in [161 [T7j . It was also shown that these Cm,n ratios are 
very stable under radiative corrections, with the LO result (including 
running of the coupling) giving very similar results to the full NLL 
calculations. 

After the arrival of LHC data it has been seen that the NLL pre¬ 
dictions, including NLO forward jet vertices, for the Cm,n ratios are in 
agreement with the experimental results. Furthermore, these observ¬ 
ables are so fine tuned to the multi-Regge limit that it is difficult for 
other approaches to fit them with accuracy. This can be seen in the 
recent studies presented in [I81II91E0], see Fig. where only a BFKL 
analysis at NLL is able to fit the large Y tail of the Mueller-Navelet 
Cm,n ratios proposed in [Ml [IT]. 

In our opinion, it is very important to continue along this line of 
work in the coming years of analysis of LHC data. In particular, it is 
needed to propose new quantities sensitive to BFKL dynamics, limit¬ 
ing the influence of the, otherwise widely dominant, collinear regions 
of phase space. In the coming Section we show that, if in the previ¬ 
ously studied events with two tagged forward jets, we tag on a third, 
central, jet there will be many distinct new observables dominated by 
BFKL dynamics which are worth investigating both theoretically and 
experimentally. 


2 Inclusive three jet production 

We propose to study events with two tagged forward jets, separated by 
a large rapidity span, and also tag on a third jet produced in the central 
region of rapidity, allowing for inclusive radiation in the remaining 
areas of the detectors. In these processes it is possible to define many 
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Figure 1: The ratio of average cosines €2,1 = C2/C1 in bins of Ay = Y, 
compared to various Monte Carlo models. Plots taken from [20] . 


differential distributions in the transverse momentum, azimuthal angle 
and rapidity of the central jet, for fixed values of the four momenta 
of the forward jets. Our predictions for these observables will use the 
BFKL formalism to describe the inclusive multi jet emission taking 
place between the three tagged jets. 

Before we proceed to explain the details of our calculation it is fair 
to highlight its limitations. We work at leading logarithmic accuracy, 
although we include running of the coupling effects. This will be im¬ 
proved in the future but we argue that, for distributions non sensitive 
to the zero conformal spin, higher order corrections will be small since 
the leading order prediction dominates the observables (the same ar¬ 
gument was used in Ref. ng I?! for the usual Mueller-Navelet jets 
case). We work with hxed four momenta for the two forward jets, this 
is for simplicity and clarity of presentation but a detailed extension 
including proper binning in these variables (which requires the intro¬ 
duction of collinear parton distributions) will be presented elsewhere 
together with further studies including the production of more than 
one jet in the central region of rapidity. Finally, although we provide 
useful analytic expressions for our observables, it is desirable to pro¬ 
duce these results using Monte Carlo event generators, a work which 
is underway. 

The two tagged forward jets A and B have transverse momentum 
kA,B, azimuthal angle 9a,b and rapidity Ya,b- The central jet is char¬ 
acterized by kj, Oj and yj and the differential cross section on these 
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variables can be written in the form 


^3^3-jet 

(Pkjdyj 


= ^ J J d^PB (pa + kj - pb) 

X p(kA,PA,YA-yj^ p(pB,kB,yj -Yb^ , 


( 1 ) 


where we assume that Ya > yj > Yb and kj lies above the experimen¬ 
tal resolution scale, p are BFKL gluon Green functions normalized to 
P (P, Q, 0) = (p- ^ and Os = asNc/n. 

In this paper we focus on quantities for which we find that the 
BFKL formalism will be both distinct from other approaches and very 
insensitive to higher order corrections. With this target in mind, we 
first integrate over the azimuthal angle of the central jet and over 
the difference in azimuthal angle between the two forward jets, Acj) = 
^A — ^B — to define a quantity similar to the usual Mueller-Navelet 

CclSG, Z.C-, 


r27T 


dA(p cos {MA(j)) / dOj 


2 -k ^3^3-jet 


d?kjdyj 


( 2 ) 



where {il) is the logarithmic derivative of Euler’s gamma function) 

P\ \ 


<t>n(p\-,P%,Y^ = 2^ ducosiuln 


PbJ njp\p% 


(3) 


Xn{B) = 2V’(1) - V’ + *2^) - V’ -(4) 


One of the experimental observables we want to highlight here corre¬ 
sponds to the mean value of the cosine of A(p in the recorded events: 


(cos (M {6a 


vr))) = 


cos(MA^)/2-d0,gf^ 


(5) 


As we have already mentioned, the perturbative stability (including 
renormalization scale dependence) of our predictions is much better 
(see m for a related discussion) if we remove the contribution from 
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the zero conformal spin which corresponds to the index n = 0 in 
Eq. (§. We can achieve this by defining the ratios 


(cos (M {6A - Ob-T^))) 
{cos {N {9 a -9b - tt))) ’ 


where we consider M, N as positive integers. 

Many studies can be performed with these ratios but here he offer 
a first one where we fix the transverse momenta of the forward jets to 
kA = 35 GeV and /cb = 38 GeV. We also fix the rapidity of the central 
jet to be one half of the rapidity difference between the two forward 
jets: yj = (Ya — Yb)I‘^ simply because this allows us to connect with 
the well-known Mueller-Navelet jets. In this way we can plot, e.g., the 
ratio TZ\ in Fig. [^for two values of the transverse momentum of the 
central jet kj = 35,40 GeV where we can see that this ratio decreases 
as a function of Ya — Ys. This is a consequence of having an increase 
in the available phase space for inclusive minijet radiation and that 
the n = 1 component decreases which energy slower that the n = 2 
contribution. 

In the BFKL formalism we have that the larger the n the slower 
the evolution with rapidity differencies. This is very important since 
it is distinct from other approaches where QGD coherence is intro¬ 
duced as it was shown in Ref. [22]. We believe this is the reason why 
the usual Monte Carlo event generators fail to properly describe the 
Mueller-Navelet ratios proposed in imiiTj and will probably also fail 
to describe the ones here investigated. 

In the last part of this work, we want to propose new observ¬ 
ables whose associated distributions have a very different behavior 
to the ones characteristic of the Mueller-Navelet case. These new 
distributions are defined using the projections on the two relative az¬ 
imuthal angles formed by each of the forward jets with the central jet, 
9a — 9j — TT and 9j — 9b — vr, in the form 



r2TT 


r2n 


d9A 


d9B / d9j cos {M {9A — 9j 


tt)) 


A3^3-jet 

COS {N {9.J -9b- tt)) ^ 

d^kjdyj 


N 

L=0 




dp^ {/) 



X 4>m {p%P^,Ya - 


\J _|_ 2y^2/j.2 pQg 

yj) (t>N {p^ + cos9,p%,yj 


(7) 
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Figure 2: A study of the ratio TZ\ as defined in Eq. (|^ for fixed values of the 
Pt of the two forward jets and two values of the pt of the tagged central jet, 
as a function of the rapidity difference between the two forward jets for the 
rapidity of the central jet chosen as yj = (Ya — Yb)/2. 


The experimentally relevant observable is the mean value in the 
selected events of the two cosines, i.e. 


(cos (M {0A - Oj - tt)) cos {N {0j -0b - tt))) (8) 

d0Ad0Bd0j cos (M {0A - 0j - vr)) cos (A {0j -0 b- 

d^Ad0Bd0j<^g^ ■ 

d^kjdyj 

As before, in order to have optimal perturbative convergence and elim¬ 
inate collinear contamination, we can remove the contributions from 
zero conformal spin by defining the ratios: 

m,n ^ (cos (M {0A -0J - tt)) cos (A {0j -0 b- tt))) 

(cos (P {0A -0J - tt)) cos {Q {0J -0b- tt))) 

and consider M, A, P, Q > 0 as integer numbers. 

We can now investigate many momenta configurations. As an ex¬ 
ample here we show some ratios 'R.^q with M, A = 1,2 fixing the 
momenta of the forward jets to kA = 40 GeV and = 50 GeV and 
their rapidities to Ya = 10 and Yb = 0. For the transverse momentum 
of the central jet we choose three values kj = 30,45, 70 GeV and we 
vary the rapidity of the central jet yj in between the two rapidities 
of the forward jets. We show the results in Fig. These distribu¬ 
tions are proving the fine structure of the QGD radiation in the high 
energy limit. They gauge the relative weights of each conformal spin 


r27r 

Jo 
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Figure 3: A study of the ratios 7^2g) 7^2g ^ defined in Eq. (|9|) for 


fixed values of the pt of the two forward jets and three values of the pt o 
tagged central jet, as a function of the rapidity of the central jet yj. 
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contribution to the total cross section. We expect the LHC data to 
agree with these results specially in the regions where yj is closer to 
(Ya — Yb)!"^- It will be very interesting to see how they compare 
to the predictions from fixed order analysis and Monte Carlo event 
generators. 


3 Summary & Outlook 

We have addressed the problem of defining new observables which can 
be sensitive to BFKL dynamics in the context of the LHC physics 
program. We have focussed on the inclusive production of three jets 
well separated in rapidity from each other. The BFKL resummation 
contains much more information than just the growth of cross sections: 
when projected on azimuthal angles it is expressed as a sum of infinite 
components from which only the first one (n = 0 Fourier component) 
grows with energy, all the others decrease. This fact has already been 
used to discriminate between this type of resummations and other 
approaches in the case of Mueller-Navelet jets, with two tagged jets 
in the very forward directions. Here we have proposed to tag a third 
jet in the central region of rapidity allowing us to test the BFKL 
formalism at a more exclusive level. This third jet connects with the 
two forward ones via two gluon Green functions and the further in 
rapidity it is emitted from them, the more decorrelated in azimuthal 
angles they will be. 

The projection on azimuthal angles we propose allows for the defi¬ 
nition of many different observables whose behavior when varying the 
pt and rapidity of the central jet is a distinct signal of BFKL dynamics. 
In order the reduce the theoretical uncertainties we have put forward 
the study of ratios of correlation functions of products of cosines of 
azimuthal angle differences among the tagged jets in Eqs. Q and (§. 
This suppresses the collinear dynamics and reduces the influence of 
the parton distribution functions. 

This program is important since there are uncertainties in the 
BFKL approach itself which need to be fixed and the current data 
recorded at the LHC will be crucial to do so. Only a fair comparison 
to experimental data can solve many of these theoretical questions. 
We believe the type of observables proposed in this work will be cru¬ 
cial to define the region of phenomenological applicability of the BFKL 
resummation. 
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